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Abstract. We study the changes in the dynamical be- 
haviour of axisymmetric spherical mean-field dynamo 
models produced by changes in their geometry and topol- 
ogy, by considering a two parameter family of models, 
ranging from a full sphere to spherical shell, torus and 
disc-like configurations, within a unified framework. We 
find that the two parameter space of the family of mod- 
els considered here separates into at least three different 
regions with distinct characteristics for the onset of dy- 
namo action. In two of these regions, the most easily ex- 
cited fields are oscillatory, in one case with dipolar sym- 
metry, and in the other with quadrupolar, whereas in the 
third region the most easily excited field is steady and 
quadrupolar. In the nonlinear regime, we find that topo- 
logical changes can alter significantly the dynamical be- 
haviour, whilst modest changes in geometry can produce 
qualitative changes, particularly for thin disc-like config- 
urations. This is of potential importance, since the exact 
shapes of astrophysical bodies, especially accretion discs 
and galaxies, are usually not precisely known. 



1. Introduction 

The magnetic fields observed in stars, accretion discs and 
galaxies are generally considered to be generated by mag- 
netohydrodynamic dynamo action. The complexity of the 
physics of such dynamos has meant that a fully self- 
consistent model is beyond the range of the computa- 
tional resources currently available, although important 
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attempts have been made to understand turbulent dy- 
namos in stars (Nordlund et al. 1992, Brandenburg et al. 
1996) and accretion discs (Brandenburg et al. 1995, Haw- 
ley et al. 1996), for example. Such studies have had to be 
restricted to the geometry of a Cartesian box, thus they 
are in essence local dynamos. However, magnetic fields in 
astrophysical objects are observed to exhibit large scale 
structure, related to the shape of the object, and thus 
can only be captured fully by global dynamo models. In 
some cases it has been possible to show that, for fully 
three-dimensional turbulence simulations, the dependence 
of the results (parity, time dependence, cycle period, etc.) 
on global properties including boundary conditions is re- 
markably similar to that of the corresponding mean-field 
dynamo models (Brandenburg 1999). This lends some sup- 
port to the validity of mean-field theory, which is used in 
the following. 

Mean field dynamo models, which attempt to capture 
such geometric features, have been employed extensively 
in order to study various aspects of the dynamics of so- 
lar, stellar and galactic dynamos (see for example, Steen- 
beck & Krause 1969, Roberts & Stix 1972). Often the 
dimensionality of the underlying partial differential equa- 
tions has been reduced by assuming axisymmetry. As the 
presence of strong differential rotation tends to destroy 
non-axisymmetric fields (e.g. Radler 1986), this restric- 
tion has some astrophysical justification, but nevertheless 
some galaxies and types of stars are known to have non- 
axisymmetric fields and so this simplification cannot be 
universally valid. 

In addition to producing solutions with pure and 
mixed parities (Brandenburg et al. 1989a, b), these studies 
have also revealed new features such as solutions with both 
chaotic behaviour (see for example Jones et al. 1985 in a 
severely truncated system, Brooke & Moss 1994, Torkels- 
son & Brandenburg 1994b, Tavakol et al. 1995) and inter- 
mittent modes of behaviour (see e.g. Brooke & Moss 1995; 
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Tobias & Weiss, 1997; Tworkowski et al. 1998, Covas et 
al. 1998a). 

An important shortcoming of these models is that they 
involve severe approximations, which essentially fall into 
two broad categories. There are those concerning the un- 
derlying physics, such as approximations involved in the 
parameterisation of the turbulent processes in solar and 
stellar convective zones (which are not known precisely), 
and there are the simplifications assumed in modelling the 
phenomenological features of these systems, such as the 
details of their geometries. These shortcomings are often 
ignored, the assumption being that changing the details 
of these models should not change their behaviours qual- 
itatively. A priori, however, there are no reasons why this 
should be so and it is therefore important that this as- 
sumption should be checked. 

Recently, an attempt was made to study the robustness 
of axisymmctric mean-field models in spheres and spher- 
ical shells with respect to changes in the functional form 
of the a effect (Tavakol et al. 1995) and also by consider- 
ing the a-quenching to be dynamic (Covas et al. 1998a). 
Here we shall employ the usual algebraic a-quenching and 
focus on the second feature discussed above, namely the 
assumptions regarding the topology and geometry of these 
models. This is important for two reasons. Firstly, these 
are the two important features differentiating the astro- 
physical situations in which such dynamos are thought to 
be operative. Thus spherical and spherical shell models are 
considered in connection with solar and stellar variability 
(with the sphere or shell representing convectively unsta- 
ble regions of stars), whereas disc (Torkelsson & Bran- 
denburg 1994a,b, 1995) and torus (Brooke & Moss 1994; 
1995) models are considered in connection with modelling 
accretion discs and galaxies. Secondly, the geometries of 
these models are often only approximations to the true 
shapes of the astrophysical systems, and so it is impor- 
tant to know the sensitivity of the dynamo solutions to 
modest changes in geometry. This is particularly true for 
dynamo models in accretion discs and galaxies. 

In this paper we study the qualitative nature of the 
dynamical behaviour of axisymmetric mean-field dynamo 
models as a function of changes in their topology, follow- 
ing transitions from sphere to shell to torus, and as a 
function of changes in their geometry resulting from alter- 
ation of the shapes of the dynamo region, whilst retaining 
the topology. We consider these changes within a unified 
framework, using the same numerical code with the same 
boundary conditions, to solve the same dynamo equations. 
In the context of disc-like models, we also consider changes 
in the rotation law. 

The models considered here include geometric features 
such as curvature, which are naturally present in astro- 
physical systems. There have also been studies of more 
idealized models in Cartesian geometry, and of low order 
models, which aim to describe generic features of axisym- 
metric dynamo solutions (e.g. Weiss, Cattaneo & Jones 



1984, Jennings et al. 1990, Jennings & Weiss 1991, Jen- 
nings 1991, Tobias et al. 1995). It is of interest to deter- 
mine whether such features persist across changes in ge- 
ometry and topology and in the highly nonlinear regimes 
considered here. We shall briefly return to this issue below 
in the light of our findings. 

In section 2 we briefly introduce the equations for the 
axisymmetric mean-field dynamo models together with 
the relevant boundary conditions. Section 3 contains a 
discussion of our results and finally section 4 contains our 
conclusions. 

2. The model 

The standard mean-field dynamo equation (cf. Krause & 
Radler 1980) is of the form 

dB 

— = V x (u x B + aB -?7 t V x B) , (1) 

where B and u are the mean magnetic field and the mean 
velocity respectively. The quantities a (the a effect) and 
the turbulent magnetic diffusivity, ?y t , appear in the pro- 
cess of the parameterisation of the second order corre- 
lations between the velocity and magnetic field fluctua- 
tions (u' and B'). For the sake of simplicity, and to fa- 
cilitate comparison with previous work, we shall ignore 
anisotropies and take a and rj t to be scalars. We assume 
that u = f2 x r = u^4>. Nondimensionalization of the 
equations in terms of a length R and a time R 2 /rjt pro- 
duces the two usual dynamo parameters C a — a R/rj t and 
Cui — floR 2 /r] t , where ao and O are typical values of a 
and |n|. For the sake of comparison with existing results 
(e.g., Covas et al. 1998a), we set |C W | = 10 4 throughout 
(except briefly, in one particular case). We use spherical 
polar coordinates (r, 9, (f) and consider axisymmetric so- 
lutions only. In order to satisfy the condition V • B = we 
solve Eq. (Q) by splitting the field into meridional and 



azimuthal components, B 



B, 



B^, and expressing 



these components in terms of scalar field functions, so that 



B, 



V x a(f>, B0 = b(f>. 



We consider the usual algebraic form of a-quenching 
namely 



ao cos 9 



(2) 



u 1 + B 2 ' 

with ao = constant. This form has been adopted in nu- 
merous studies. 

In our calculations we consider two types of rotation 
laws, one with constant shear and the other with approxi- 
mately rigid rotation close to the rotation axis and tending 
to Keplerian rotation at larger radii. Hereafter we refer to 
the latter as a "generalized Keplerian" rotation law. The 
first constant shear profile has been used extensively in 
studies in many different geometries. We recognise that 
helioseismological evidence indicates that the rotational 
profile of the sun gives a shear in a very concentrated re- 
gion. Our main aim is to allow comparison with a range of 
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published models, hence we chose this simple and widely 
used form. The second profile has been extensively used 
in studies of accretion discs. In that it gives a differential 
rotation that decreases with radius it might be considered 
to give a generic representation of the rotation in a galac- 
tic disc, although rotation profiles in galaxies are certainly 
not Keplerian. 

For the rotation profiles that we will consider here, we 
take 



u = Qvcrtp, 



(3) 



where w = rsin# is the distance from the rotation axis. 
For the first profile mentioned above, we shall, consistent 
with Brandenburg et al. (1989a), take Q — £lo r , where fio 
is constant, and for the second profile, we shall take 



n = n 



1 



w 
— o 



-1/r, 



(4) 



where q — 3/2 for accretion discs. In conformity with 
Torkelsson & Brandenburg (1994a), we take vjq to be 0.3 
and n = 10. Relation approximates solid body rota- 
tion for small w and approaches Keplerian rotation for 
zu ^> -cuq. Note that for this profile dVl/dr < when 
C w > 0, whilst for profile (||) with constant shear the 
signs are reversed, i.e. dVl/dr > when C u > 0. 

The models we shall consider here are all constructed 
from a complete sphere of radius R. We can remove an in- 
ner concentric sphere of radius ro (in order to go smoothly 
from a sphere to a spherical shell). The change from a 
sphere or a spherical shell to a torus- or disc-like config- 
uration is produced by excising a cone of semi-angle 9q 
about the rotation axis, from the north and south polar 
regions. The configuration is illustrated in Fig. [I], where 
the dynamo region is produced by revolving the merid- 
ian section shown about the vertical dotted rotation axis 
through the origin O. 

In order to solve the dynamo equations, boundary con- 
ditions need to be imposed on the surfaces of the dynamo 
region created by removing the cones (denoted by S c ), on 
the inner surface r = tq (denoted by Si), and also on the 
outer curved boundary (denoted by S ), as depicted in 
Fig. | 

We adopt a perfect conductor boundary condition, 



0, 



d(rb) d(ra) 



dr 



di 



Si 



(5) 



Other choices are possible (see, e.g., Tworkowski et al. 
1998), but this has been used commonly in earlier inves- 
tigations. 

On S a vacuum boundary condition has often been 
used. However this, although appealing, is actually quite 
arbitrary as, for example, there is no evidence that the 
solar poloidal field is curl-free immediately above the sur- 
face. When 9 > 0, implementing an essentially non- 
local vacuum boundary condition is not straightforward, 




O 



Fig. 1. A schematic representation of the dynamo region 
produced by the cuts characterised by the model param- 
eters ro and 9q. 



and so we have chosen a plausible local condition. Other 
authors have previously taken a radial field condition, 
d{ar)/dr = 0, (e.g. Gilman & Miller 1981). We used 



da/dr = 0, b = 0, 



(6) 



which results in a poloidal field on S which is typically 
near to radial over most of the boundary. We stress that, 
in the absence of a detailed model of processes above a 
stellar surface, all choices of boundary condition on S 
are more-or-less arbitrary, within general restrictions of 
plausibility. 

If a perfect conductor boundary condition is used on 
the surfaces S c then as 9o — *• the model does not tend 
to that used for the full sphere, where a = b = on the 
axis. Thus we set 



on S c , 



(7) 



to ensure continuity of behaviour as 0q — > 0; again other 
choices are possible. 

In the following, as is customary, we shall discuss the 
behaviour of the solutions by monitoring the total mag- 
netic energy, E = tJ^ J B 2 dV, taken over the dynamo 
region given by r < r < R and 9o < < n — 9 . We 
split E into two parts, E = E^ + E (s \ where 
and E( s > are respectively the energies of those parts of 
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the field whose toroidal field is antisymmetric and sym- 
metric about the equator. The overall parity P is given 
by P = [E^ - E^]/E (Brandenburg et al. 1989a), so 
P = — 1 denotes an antisymmetric (dipole-like) pure par- 
ity solution and P = +1 a symmetric (quadrupolc-likc) 
pure parity solution. 

For the numerical results reported in the following sec- 
tions, we used a modified version of the axisymmetric dy- 
namo code of Brandenburg et al. (1989a). We took a grid 
size of 41 x 81 mesh points in the dynamo region. To test 
the robustness of the code we verified that no qualitative 
changes were produced by employing a finer grid, different 
temporal step length (we used a step length of 10~ 4 R 2 /rj t 
in the results presented in this paper). We considered a 
family of models ranging from a full sphere through spher- 
ical shells to torus- and disc-like configurations. 

3. Results 

We present the results of our comparative studies of dy- 
namo solutions in regions that are produced by variations 
of ro and do m the ranges < tq < 1 (in dimensionless 
units) and < 8q < ir/2 respectively. In particular, we 
took the following values: 

r = 0, 0.2, 0.5, 0.7 (8) 
O = 0, 4.5°, 22.5°, 45°, 75° (9) 

To present our results, we use the following notation in 
the figures: the prefix "steady" denotes a constant energy, 
whilst "A" and "S" on their own respectively represent 
pure antisymmetric (dipolar) (P = —1) and symmetric 
(quadrupolar) (P = +1) solutions with periodically os- 
cillating energy, "OM" represents solutions which possess 
periodic oscillations in both parity and energy, and "C" 
and "I" denote chaotic and intermittent behaviours re- 
spectively. By intermittent behaviour we loosely mean ir- 
regular switchings between two or more different dynam- 
ical regimes characterised by different statistics. (For a 
precise characterisation of this type of behaviour in some 
dynamo models see Covas et al. 1998a.) We should add 
here that the computational cost of our calculations lim- 
ited the resolution to which we were able to search the 
parameter space. Consequently all our results regarding 
the sequence of parity changes are subject to this limita- 
tion. 

Since many previous authors have discussed behaviour 
at the onset of dynamo action, we shall present our results 
in two subsections. First we present our results concern- 
ing the onset of dynamo action as we go from a sphere 
through spherical shells to torus- and disc-like structures, 
and then we consider the dynamical behaviour for each 
configuration as the dynamo number is increased. Our de- 
tailed results are presented in Table and Figs. | to [If 



3.1. The onset of dynamo action 

We summarise our onset results by considering the two 
types of deformation of the basic sphere in turn. Starting 
with the case where an inner sphere is removed, giving 
a spherical shell, we note that Roberts (1972) found nu- 
merically (but for different boundary conditions) that for 
an aid dynamo with a negative dynamo number, the most 
easily excited mode has dipole-like symmetry for < 0.7, 
whereas for ro > 0.7 the most easily excited mode has 
quadrupole-like symmetry. Our corresponding results are 
summarised in Fig. |2[ In particular, the top row of this fig- 
ure, summarizing the behaviour of a sphere and spherical 
shells without polar cones removed, agrees qualitatively 
with the results of Roberts (1972). 



r =0.0 r =0.2 r =0.5 r =0.7 




Fig. 2. Schematic representation of the dynamo regions 
and their corresponding onset parity with = —10 and 
a constant rotational shear profile. 



The next three rows in this figure correspond to the re- 
moval of cones of semi-angles 8q = 4.5°, 22.5° and 45° re- 
spectively. These results indicate that despite the removal 
of the polar cones, the onset behaviour still remains the 
same as that given in Roberts (1972). In this way we may 
say that the onset results for spheres are quite robust with 
respect to the polar cuts considered here. 



Covas et at: Mean- field dynamos with different geometries and topologies 



5 



Steady S Steady S Steady S 



10 





Steady A Steady A Steady A Steady A jc^+io 4 
Steady A Steady A Steady A Steady A Jc^-io 4 



in 
o 
II 

CD 





Steady S Steady S Steady S Steady S Jc^+io 4 

Fig. 3. Schematic representation of the dynamo regions 
and their corresponding onset parity. The top and bottom 
panels are for constant shear and generalized Keplerian 
rotational profiles respectively. 



Table 1. Values of C a for the onset of dynamo action 
as a function of the shell radius ro and the cut angle 6*o 
with = — 10 4 . The row labeled * is produced using 
the generalized Keplerian rotation profile given by Eq. (^J) 
with C u = — 10 4 . Here SS indicates a Steady S state. 






r = 0.0 


r = 0.2 


r = 0.5 


r = 0.7 





0.45 A 


0.38 A 


0.31 A 


0.42 S 


4.5° 


0.45 A 


0.38 A 


0.32 A 


0.42 S 


22.5° 


0.70 A 


0.60 A 


0.42 A 


0.50 S 


45° 


1.7 A 


1.3 A 


0.9 A 


1.0 S 


48° 


1.86 A 


1.48 A 


1.04 A 


1.07 S 


49° 


1.98 A 


1.58 A 


1.10 A 


1.113 A 


50° 


2.1 A 


1.679 A 


1.15 A 


1.17 A 


55° 


3.0 A 


2.37 A 


1.56 A 


1.493 A 


60° 


4.45 A 


3.545 A 


2.27 A 


1.98 A 


65° 


7.2 A 


5.72 A 


3.57 A 


2.77 A 


67° 


9.5 A 


7.1 A 


4.4 A 


3.24 A 


70° 


11 SS 


8.05 SS 


3.09 SS 


4.295 A 


75° 


13.0 SS 


10.1 SS 


5.8 SS 


7.6 A 


75°* 


4.0 SA 


3.7 SA 


3.5 SA 


11.5 SA 



To compare our results further with those from in pre- 
vious studies, for example of dynamo action in tori by 
Brooke & Moss (1994, 1995) and in discs by Torkelsson & 
Brandenburg (1994a, b, 1995), we also studied the case of a 
large polar cut angle (8q = 75°), with both constant shear 
and generalized Keplerian rotation profiles, and these re- 
sults are summarized in Fig. |3|. Note that, as commented 
after Eq. (@), with the generalized Keplerian rotation law 
(0), C w > corresponds loosely to the constant shear law 



with C u < 0. However equal values of \C^\ do not imply 
equal values of the rotational shear. 

As can be seen from a comparison of Figs. || and |[ 
the onset behaviour in these models changes from oscil- 
latory to steady solutions as the cut angle 6>o increases. 
Such a change in behaviour has also been found in oblate 
spheroidal geometry as a sphere flattens to a disc-like 
configuration (Stix 1975, Soward 1992a, b and Walker & 
Barenghi 1994). According to these authors, the crucial 
factor in the change from oscillatory to steady behaviour 
is the aspect ratio of the spheroids, defined as the ratio of 
the minor to the major axis. They speculated that this was 
because dynamo waves were restricted from propagating 
in the vertical direction (i.e. parallel to the minor axis). In 
this context it is interesting to note that in our case the 
only oscillating solution for 6*o = 75° occurs when ro = 0.7 
and the above aspect ratio is O(l), as for a sphere. 

In Table | we summarise the values of C a correspond- 
ing to the onset of dynamo action in each configuration. 
Except for the case with generalized Keplerian rotation 
law, this shows that for a fixed ro the values of C a at 
the onset of dynamo action increase as the cut angle 6*o 
increases for almost all configurations as the shells get 
thicker. 

We should also note that the steep increase in the onset 
value of C Q , as we go from 9q = 45° to Qq — 75°, may 
be related to the transition from shell-like to disc-like 
configurations. Table [j] in conjunction with Figs. || and 
H confirm that in the constant shear case configurations 
with smaller polar cut angles follow the linear results due 
to Roberts (1972), while those with larger cut angles (9o i<J 
49°) do not, thus suggesting the possible existence of a 
transition as the cut angle is increased. 

For the disc-like configurations (with 6>o ib 70°), the 
onset of dynamo action is to a steady state when ro < 0.7, 
when the symmetry is determined by the sign of If 
this is positive, that is for the constant shear with C w = 
+10 and for the generalized Keplerian rotation profile 
with C u = — 10 4 , the onset modes are antisymmetric. For 
the negative case the onset mode is always symmetric. 

3.2. Nonlinear regime 

We present our results in the nonlinear regime in three 
stages. First we consider the transition from the full sphere 
to a spherical shell, then from a full sphere to that with 
axial cones removed and finally we study a full sphere 
with both an inner core and axial cones removed simul- 
taneously. Our detailed results are presented in Figs. [I] - 
||, corresponding to five different values of 8q. Each con- 
tains four panels for each of the values of ro considered, 
the panels from top to bottom corresponding to ro val- 
ues of 0.7,0.5,0.2 and 0.0 respectively. We have concen- 
trated mostly on the more usual case of negative dynamo 
number (C w = — 10 4 ), but for comparison we also con- 
sider, for large 9 , the case of a positive dynamo number 
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(C w = +10 4 ; see Fig. ^|), as well as the case with general- 
ized Keplcrian differential rotation, which is relevant for 
accretion discs (Figs. ITol-O). 
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Fig. 6. Partial bifurcation diagram, as Fig. [|, for the case 
9 = 22.5° with constant (negative) shear, C w = — 10 4 . 



Fig. 4. Partial bifurcation diagram, showing the energy in 
the stable solution, for the case of complete spherical shells 
(#o = 0), with constant (negative) shear, C u = — 10 4 . 
For oscillatory modes, the mean energy is plotted. Here 
and in following figures, the panels from top to bottom 
correspond to r values 0.7,0.5,0.2 and 0.0 respectively. 
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Fig. 5. Partial bifurcation diagram, as Fig. [|, for the case 
of a polar cut with a very small angle 9q = 4.5° and with 
constant (negative) shear, C w = — 10 4 . 



Before discussing our results in detail we emphasise 
that the sequence of parity changes reported below are 
those that we found to within the resolution of the C a pa- 
rameter employed here. It is therefore important to bear 
in mind that unless there is a clear theoretical reason for 
believing in the presence of a particular sequence of par- 
ity changes in the nonlinear regimes, it is always possible 
to miss finer details of such sequences for any given fi- 
nite resolution in C a . Our standard step in C a is 0.5, but 
on occasion this was very much reduced (see, e.g., Fig. 
|l3| ). In addition, although we have attempted to integrate 
over times long enough to eliminate any transient effects, 
nonetheless, it is not possible to always rule out the pres- 
ence of very long-lived transients. With these qualifica- 
tions we proceed to discuss our results in more detail. 

The results concerning the transition from a full sphere 
to spherical shells are shown in Fig. [|. This shows that 
the transitions, from the onset of dynamo action into the 
nonlinear regimes as C a increases, are the same for the full 
sphere and the two thickest shells (i.e. ro = 0.0,0.2,0.5). 
All possess the same initial sequence of parity changes 
given by 



A, OM, S, A. 



(10) 



The only difference is that, as the shells become thinner 
(ro increases), the sequence is compressed to a narrower 
interval of C a . For the thinnest shell (tq — 0.7), however, 
the initial sequence is 



>S\ A, S, A. 



(11) 



When a polar cone is removed (6q > 0), the initial 
sequence of parity changes is the same for the full sphere 
and the two thickest shells for 6q < 45°. The behaviour for 
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the thinnest shell (7-0 = 0.7) and the largest polar cut angle 
(8q = 75°) differ. We also observe that the removal of the 
polar cones causes more changeable and varied behaviour, 
as does the removal of the inner cores, providing that the 
first bifurcation is to an oscillating solution. 
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Fig. 7. Partial bifurcation diagram, as Fig. ||, for the case 
8q = 45°, with constant (negative) shear, C u = — 10 4 . 

Figures || to (?] show that there are similarities in the 
initial sequences of parity changes as 80 changes. The con- 
figurations with 6» = 4.5°, 22.5° and 45° and r = 0.0, 0.2 
and 0.5 all have the same initial sequence, namely, 



A, OM, S, A, 



(12) 



whilst configurations with tq = 0.7 and these same angles 
have the sequence 



S, A , A. 



(13) 



These two initial sequences agree with the corresponding 
ones in the full sphere case (8q = 0°). Differences in the 
bifurcation sequences tend to appear at larger values of 
C a where the dynamo becomes more non-linear. Also at 
these higher values of C a , more exotic behaviours such 
as intermittency and chaos can appear, for example for 
O = 0° and 4.5° with r = 0.2 and 0.5, 8 = 22.5° with 
tq = 0.2, 0.5 and 0.7, and 80 — 45° for all values of ro 
considered here, whilst chaotic and intermittent behaviour 
only occurs for 8q = 75° when ro = 0.7. 

When O = 75° (Fig. g) the behaviour is totally differ- 
ent. The steady quadrupole mode (S) remains stable for 
a very large range of values of C a . At some point an oscil- 
latory S mode takes over. For r = 0.5 OM modes appear 
and only for ro = 0.7 are oscillatory modes observed for 
all values of C a studied. 
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Fig. 8. Partial bifurcation diagram, as Fig. |, for the 
case of a large polar cut angle, O = 75°, with constant 
(negative) shear, C w = — 10 4 . 
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Fig. 9. Partial bifurcation diagram, as Fig. f§ for the 
case of a large polar cut angle, 8q — 75°, with constant 
(positive) shear, C w = +10 4 . 



We also studied the effects of changing the sign of 
on our results and this is shown in Fig. ^|. As can be seen, 
a steady A mode is first excited for all r . The general 
behaviour is also less complex and more uniform, partic- 
ularly if we compare the panel corresponding to ro = 0.7 
for both signs of C w . 

The effects of changing the rotation law to the gener- 
alized Keplerian rotation law are depicted in Fig. |lC| and, 
as can be seen, the steady S mode and the subsequent 
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oscillatory S mode are stable for all values of C a consid- 
ered. Only in the case ro = 0.7 is new behaviour seen, in 
that there is an interval (C a as 15) where only the triv- 
ial solution with zero magnetic field exists. Such a gap in 
the bifurcation diagram has been found earlier in connec- 
tion with torus and accretion disc dynamos (e.g. Brooke 
& Moss 1995, Torkelsson & Brandenburg 1994b, see also 
Ruzmaikin et al. 1980 and Stepinski & Levy 1991 for lin- 
ear results). Immediately after this gap an antisymmetric 
mode appears. The existence of this gap does not appear 
to be a peculiarity associated with the particular value of 
we used, as we also found it when C w = 2.7 x 10 4 , a 
value chosen to give the same value of the rotational shear 
at the mean radius of the shell (r = 0.85) as the constant 
shear model with C w = — 10 4 , shown in Fig. @. We show 
these results in Fig. [ll]. It is clear from comparing Fig. [io] 
and Fig. [ll] that the bifurcation sequence is sensitive to 
the value of C u , but the gap in which no field is excited 
remains. 



10 

A 10' 

I 10 ° 
V lir 1 

10" 2 
10 2 

A 10' i 

I io" t 
v io"' 

10" 2 

io 2 
A 10 ' 

v 10"' 

10" 2 

io 2 

A 10' 

v io- ] 
10 : 



Steady S 



Gap 



Steady S 



SlL.nl> 





Steady S 








S j 



■ Steady S 

• Oscillating A J 
' Oscillating S 

• Oscillating M 

■ Chaotic 



50 
C„ 



Fig. 10. Partial bifurcation diagram, as Fig. ||[ for the 

case of a large polar cut angle, 9q = 75°, with the gener- 
alized Keplerian rotation law, C u = +10 4 {dVt/dr < 0), 
= 0.3. 



When the sign of the generalized Keplerian rotational 
shear is reversed (Fig. |l2|), a steady A mode is again first 
excited, similar to the case of linear shear. Only for ro = 
0.7 is the behaviour slightly different in that there is a S 
mode that appears for C a as 22. 

In the above figures, the nature of the solutions (and 
thus the symbols shown) are all determined only to within 
the numerical resolutions chosen and for the length of of 
time during which solutions were followed and thus in 
which transients could decay. In particular, near the bi- 
furcation points (as in the case of transitions from pure 
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Fig. 11. Partial bifurcation diagram, as Fig. ^, for the 
case of a large polar cut angle, 0$ = 75°, with the general- 
ized Keplerian rotation law, C u = +2.7xl0 4 {d^l/dr < 0), 
zuq = 0.3. 
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Fig. 12. Partial bifurcation diagram, as Fig. ^, for the 
case of a large polar cut angle, 9 = 75°, with negative 
generalized Keplerian rotation law, C u = — 10 4 (dfl/dr > 
0), uj = 0.3. 



parity A to S (or vice versa)), we find as expected ex- 
tremely long transients (at times of more than hundreds of 
diffusion times). Bearing all this in mind, we do find some 
evidence which suggests that the conjecture by Jennings 
& Weiss (1991), according to which transitions from one 
pure parity to the opposite are mediated by a mixed parity 
regime, may not always hold. An example of this type of 
transition is given in Fig. [l3], where we have used an excep- 
tionally fine resolution in C a . In this case we found after 
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some hundreds of diffusion times that with C a — 4.2868 
and C a = 4.2869 the parity became P = +0.9999633 and 
P = —0.9999967 respectively. A conclusive resolution of 
this question requires a more substantial study, and we 
intend to return to this point in a future publication. 
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Fig. 13. Parity versus time diagram for two runs with 
same initial conditions except the value of C a . The other 
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Our results also indicate that overall the behaviour 
appears to be more complicated in shells than in the full 
sphere dynamo models, a feature also observed in Covas 
et al. 1998a. 

Also, the nonlinear solutions in configurations which 
correspond most closely to a thick torus, i.e. 9 = 22.5°, 
r = 0.2, 0.5 and also 6 — 45°, r — 0.5, display chaotic 
behaviour at the lowest values of C a . This corresponds 
with the findings of Brooke & Moss (1994, 1995), that 
chaotic behaviour was excited at modestly supercritical 
dynamo numbers in tori, whereas chaotic behaviour only 
occurred at highly supercritical dynamo numbers in disc 
and shell configurations. 

To obtain an overview of the behaviour of the two pa- 
rameter family of models considered here, we recall that 
Table [j] suggests the presence of three different regions 
with distinct onset behaviours in the parameter space 
(ro,6*o). In the constant shear rotation law case two of 
these regions possess oscillatory dipolar and quadrupolar 
onset symmetries respectively and in the third the most 
readily excited field is steady and quadrupolar. These re- 
gions could also be interpreted from a different point of 
view by considering the behaviour observed in previous 
dynamo models. Viewed in this way, there are at least 
two transitions in the (ro,#o) parameter space. The first 
corresponds to a transition from a region which is consis- 



tent with Roberts' (1972) results (characterised by an A 
onset mode for ro < 0.7 and by an S onset mode with 
fo > 0.7) to one that does not (i.e. one with antisymmet- 
ric onset modes for all ro). The second transition is from 
a disc-like regime in which steady modes are commonly 
found to one in which they are absent. In order to eluci- 
date this second transition we have made two studies of 
the parameter space at a finer resolution, by fixing ro and 
varying 9q and vice versa. These results are depicted in 
Figs. |TJ and [l5| respectively at fixed values of C a . (We 
note that changing the value of C a does not qualitatively 
change our conclusions.) Those results indicate the pres- 
ence of two sharp transitions in mean energy and/or parity 
at around 9q ~ 70° for r ~ 0.5 and O ~ 75° for r ~ 0.6. 
These features are also corroborated by the time series of 
energy and parity. 




Fig. 14. Transition from sphere/shell to disc-like be- 
haviour for a constant shear profile as a function of 9q 
for C a = 10, r = 0.5, C w = — 10 4 . The inset gives the 
corresponding behaviour of the parity. 

Finally, we have also examined the behaviour of our so- 
lutions as a function of the aspect ratio, which we defined 

as 



.4 



2(1 -r ) 



(l + ro)(7r-20 o )' 



(14) 



With this definition, the smallest aspect ratio corre- 
sponds to the configuration with r = 0.7 and 6>o = 
(A = 0.11), and the largest to that with ro = and 
6> = 75° (A = 3.82). We find that at large values of 
C a , for example C a > 30, the averaged energy (E(t)) per 
unit volume of the dynamo region is largest for the con- 
figuration with the smallest aspect ratio, and vice versa. 
In addition, for these values of C a , we find that the aver- 
age energy per unit volume decreases as 9q increases (i.e. 
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Fig. 15. Transition from disc-like to planar like behaviour 
for a constant shear profile as a function of ro for C a — 60, 
8q = 75°, C w — — 10 4 . The inset gives the corresponding 
behaviour of the parity. 

the cuts get larger) and as r decreases (i.e. the shells get 
thicker). 

4. Conclusions 

We have made a comparative study of the dynamical be- 
haviour of a two parameter family of axisymmetric mean- 
field dynamo models, ranging from a full sphere to shells 
to torus- and disc-like configurations, within a unified nu- 
merical framework. We have found overall agreement with 
previous studies of axisymmetric mean-field dynamos for 
spheres, shells, tori and discs. More precisely we find the 
following: 

1. Our results show that changes in both topology and ge- 
ometry affect the dynamical properties of our models, 
although these two factors arc inevitably interlinked. 
Specifically, we have three topologically distinct con- 
figurations, sphere-like, shell-like and torus-like. Each 
of these topological configurations can undergo a con- 
tinuous sequence of geometric changes, from flat ob- 
jects where the vertical extent is much less than the 
horizontal, to objects where the vertical and horizon- 
tal extents are similar. Our principal findings are that 
torus-like topologies seem to produce more complex 
behaviour and that, within a given topological configu- 
ration, modest changes in geometry can give rise to dif- 
ferent bifurcation sequences. These findings may have 
considerable astrophysical significance, in that the ac- 
tual shapes of galaxies and accretion discs often do not 
conform to the simple geometric shapes used in many 
studies of mean-field dynamos. Also, since our study 
concerns the geometry and topology of the dynamo- 
active region, restrictions of the dynamo-active region 



in stellar convection zones may have important con- 
sequences for the behaviour of the dynamo and these 
should be taken into consideration when comparing the 
predictions of dynamo models with observed magnetic 
behaviour in astrophysical objects. 

2. We confirm that the full sphere and shell results of 
Roberts (1972), concerning the symmetry of the on- 
set modes, are robust with respect to changes in the 
outer boundary condition, and we extend them to the 
configurations with very small to medium polar cut an- 
gles. Steady modes are found for the largest polar cut 
angles, agreeing with results from oblate spheroidal ge- 
ometry. These modes are robust under a change in the 
rotation law, as found by Walker & Barcnghi (1994), 
except in the case where the configuration resembles 
a thin torus rather than a disc. In this case an oscil- 
lating dipolar solution occurs at onset for a constant 
shear rotation law, but a steady quadrupolar solution 
is found with a generalized Keplerian law. In the latter 
case the steady solution disappears at higher values of 
C a and there is a gap before an oscillating antisym- 
metric solution appears. Such a gap was previously 
observed with a generalized Keplerian rotational pro- 
file in both disc (Torkelsson & Brandenburg 1994b) 
and torus geometry (Brooke & Moss 1995). Stepin- 
ski & Levy (1991) also found a similar gap in linear 
theory. It should be noted that the presence of a gap 
does depend on the nature of the nonlinearity: it is 
absent if magnetic buoyancy is the dominant nonlin- 
earity (Torkelsson & Brandenburg 1994a). Our results 
here indicate the steady solution is also fragile with 
respect to changes in the rotation law. 

3. For the largest polar cut angle, we find evidence for 
a different mode of behaviour in the thinnest shell 
(ro = 0.7) from that found in the thicker shells and 
full sphere, namely a more complicated behaviour with 
denser bifurcation sequences, including chaotic and in- 
termittent solutions. Apart from the thinnest shells 
(ro = 0.7) the dynamics is qualitatively the same for 
each sign of the shear. This is consistent with the ob- 
servation that the thinnest shell approximates a setting 
in which curvature plays a smaller role. 

4. We note a consistent finding in our results, that 
dynamo-active regions corresponding topologically 
and geometrically to a torus exhibit chaotic behaviour 
at lower dynamo numbers compared to spherical, shell- 
like or disc-like configurations. This concurs with the 
results of Brooke & Moss (1994,1995), who found 
transition to chaotic behaviour for a dynamo in a 
toroidal volume at dynamo numbers that were ap- 
proximately three times supercritical. Given that the 
solar (and possibly some stellar) magnetic cycle ex- 
hibits irregularity on a variety of time scales, and that 
many authors have speculated that such behaviour is 
a manifestation of deterministic chaos, this result may 
have some astrophysical significance. It therefore seems 
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worthwhile to investigate such dynamos, with a variety 
of rotational profiles, possibly corresponding to differ- 
ent types of stars, and using different forms of nonlin- 
earity to saturate the field. Since it is the restriction of 
the dynamo active volume to a toroidal configuration 
that appears to be important, these dynamos could be 
considered as embedded in convective shells represent- 
ing the outer regions of stars. 
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